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ABSTRACT 


The  first  part  of  the  report  contains  a general  parametric  analysis  of 
the  heeling  angle  of  a submarine  exposed  to  currents.  The  current  is  assumed 
to  have  an  arbitrary  velocity  distribution  with  respect  to  depth  and  to  lay 
within  a plane  perpendicular  to  the  length  axis  of  the  submarine. 

In  the  second  part,  specific  observed  current  conditions  from  Hood  Canal, 
Washington  are  applied.  Records  from  a total  of  ten  current  meters  are 
evaluated  and  applied  in  computing  cumulative  histograms  of  the  heeling  angle. 
The  "it"  and  "v"  co-ordinate  system  for  the  velocity  components  is  rotated 
(relative  to  the  true  North)  by  two  angles  (-10  and  +50  degrees  angular)  to 
reflect  actual  location  of  the  mooring  piers.  Also,  the  effect  of  both  flood- 
ing and  ebbing  currents  is  evaluated  by  presenting  the  analysis  for  both  posi- 
tive and  negative  velocity  components. 
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J.G.  Dworski  : 

HEELING  ANGLE  OF  SUBMARINE  MOORED  NORMAL  TO  CURRENT  FLOW 


Introduction 

"Trident"  class  submarines  will  be  moored  at  open  piling  piers  in  Hood 
Canal  and  will  - therefore  - be  exposed  to  fairly  strong  currents.  A project 
contractor's  study  proposes  a procedure  whereby  the  current  velocity  would 
be  measured  at  three  depths  (5  m,  10  m,  15  m),  a quadratic  function  fitted 
to  these  data,  and  the  heeling  moment  computed  considering  the  submarine 
hull  as  a vertical  flat  plate  of  same  area  as  the  vertical  projected  area 
of  the  actual  hull. 

Such  procedure  will  overestimate  the  heeling  moment  because:  as  the 
actual  cylindrical  surface  of  the  submarine's  hull  curves  away  from  the 
vertical,  there  is  - intuitively  - less  contribution  by  the  "normal  drag" 
and  more  contribution  by  the  "skin  friction"  drag.  The  latter  is,  however, 
comparatively  small. 

This  research  note  decomposes  the  drag  forces  into  a component  normal 
to  the  hull  surface  and  a component  tangential  to  it.  Also,  it  uses  some 
experimental  data  - curve  fits  - developed  by  the  U.S.  Naval  Ordinance  Test 
Station  in  China  Lake,  California,  for  drag  on  inclined  surfaces,  referred 
to  experiments  determining  drag  in  normal  flow  incidence.  The  resulting 
heeling  moment  is  then  computed  in  general  vectorial  form  and  covers  any 
placement  of  the  restraining  cleats,  linking  the  submarine  to  the  pier.  Al- 
so, different  levels  of  displacement  (immersion)  of  the  submarine  are  para- 
metrically accounted  for. 

The  main  result  is  that  the  actual  heeling  angle  will  be  only  about  2/3 
of  that  predicted  by  the  "flat  plate"  approximation.  A "Trident"  class  sub- 
marine will  heel  over  about  5°-6°  in  2 kt  of  current  if  restrained  at  cleats 
flush  with  the  submarine's  deck. 


Analysis 

Geometrical  Constraints 

Geometrical  constraints  are  shown  in  Fig.  1. 
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a = cos 

Vector  from  cleat  to  £ hull,  C ^ (or  C or  Cfl ) is  fixed  vector, 
for  example 

C = (d  + p - R)k  = Dfe 
o 1 a o 

because  submarine  heel  angle  is  assumed  to  be  small. 


Normal  and  Tangential  Forces  for  Decomposition  Model 

The  following  model  has  been  developed  in  a study  made  for  Honeywell 
Defense  Laboratory,  Seattle,  by  Dworski  (1965).  It  is  based  on  data  collect- 
ed by  the  U.S.  Naval  Ordinance  Test  Station,  China  Lake,  California  (1962) 
relating  drag  on  inclined  surfaces  to  wind  tunnel  and  flume  tests  determin- 
ing drag  under  normal  (perpendicular)  incidence  of  flow  upon  the  surface. 

The  empirical  formula  determined  therefrom  was  modified  in  the  above  cited 
report  to  include  certain  compatibility  relations  that  become  important 
at  large  inclinations  (where  accurate  description  of  tangential  forces  be- 
comes mandatory).  The  appropriate  schematic  is  shown  in  Fig.  2. 


Figure  2 


Let  G - Cp  ^ U\ll\  (force  per  unit  area)  (6) 

where  C ^ is  the  "normal"  drag  coefficient,  p the  density, 
and  where  |i/|  is  the  absolute  value  of  U which  is  to  be  used  so  that  U(z) 
may  be  defined  negative  for  some  z (i.e.,  the  flow  may  reverse  at  some 
depth) . 


Dworski,  J.G.  (1965)  "Mathematical  Analysis  of  Tow  Cable  Configuration" 

Honeywell  Defense  Laboratory,  Seattle 

Project  "Nixie"  and  Project  "Lolita",  Classified  (above  segment  detached  from 
Project  report  and  declassified) 
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where  n - 0.02  is  the  friction  coefficient  on  a flat  plate  (a 
function  of  the  Reynolds  number  and  roughness)  multiplied  by  C^. 

Note  that,  in  Eq  (7)  and  (8),  the  absolute  values  |s-£n0|  , |cos0|  per- 
mit the  angle  of  inclination  0 to  go  beyond  tt/2. 

Consider,  for  better  visualization,  the  limiting  cases: 

from  Eq  (7)  and  (8):  Table  1 
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At  any  intermediate  angle  0 < 0 < ir/2  or  tt/2  < 0 < ir  the  formula- 
tion of  Eq  (7)  and  (8)^_affords  the  proper  analytical  description  of  the 
vertical  component  of  F ' and  F ^ as  controlled  by  the  sign  of  the  product 
sin  0 cos  0 ; the  horizontal  components  are  — because  of  the  squared  terms 
sin^Q  and  cos20  — positive  as  long  as  U is  positive.  Eq  (6)  permits  - 
as  stated  before  - that,  at  some  z,  U becomes  negative  affording  proper 
description  of  any  arbitrary  horizontal  flow  field. 


General  Case 

Computation  of  the  "heeling"  moment  (of  flow  forces  on  the  submarine 
hull)  relative  to  any  fixed  restraining  point  (cleat  point)  is  facilitated 
by  the  vectorial  decomposition  shown  in  Figure  1,  decomposing  that  moment 
into  two  integrals  for  each  of  the  force  subsets  ( F ' and  ?y) : 

(a)  moment  relative  to  the  center  of  the  cylindrical  hull,  which 
is  zero  for  the  normal  forces  . 

(b)  moment  relative  to  the  fixed  point  - related  to  the  center  of 
the  hull: 

Consider  a unit  width^Aj/  = 1 ft)  strip  of  the  submarine  hull  (the  y 
axis,  and  the  unit  vector  j associated  with  it,  aia  positive  out  of  the 
plane  in  Fig.  1;  moments  are  positive  when  inducing  a mathematically  positive 
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(counterclockwise)  rotation  about  the  fixed  point) 


Then  the  heeling  moment  relative  to  Cj  or  Cg  (or  C ) is  given  by  the  vector 
(outer  or  cross)  product  integrals: 
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Note  that  the  current  shear  is  reflected  in  description  of  and 

through  Eq  (6)  in  conjunction  with  Eq  (2).  For  example:  if  the  current 

profile  were  given  as  a quadratic  function  of  depth: 
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U(z)  = aQ  + a.jZ  + a^z  (10) 

Then  using  Eq  (2),  within  reach  of  the  submarine  draft: 

U(8)  = aQ  + dj  [d  - R(1  - cos  Q)  + a^d  - R(1  - cos  0)]2  (11) 
so  that  for  6 = 0,  U(Q=0)  = U(d)  = aQ  + a^d  V-  a^d2  (12) 

and  for  0 = ir-a,  U(d= u-a)  = U(o)  = a.Q  (13) 

where  Eq  (13)  follows  from  the  fact  that  by  Eq  (4):  cos(-n-a)  = 1 - d/R 

In  general,  whenever  the  current  has  a z dependence,  Eq  (9)  is  best  evalu- 
ated by  a small  computer  program  to  avoid  the  messy  algebra. 

If  the  current  were  independent  of  z one  can  easily  evaluate  Eq  (9)  in 
analytic  form,  as  shown  in  the  following  specialization. 
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Special  Case,  Analytical  Solution 

Assume  current  to  be  independent  of  depth  and  assume  the  cleat  position 
to  be  at  the  center-plane  of  the  submarine  (i.e.,  at  C ) . Evaluate  the 
heeling  moment  Mr  term-by-term,  using  R(Q)  as  defined  Sy  Eq  (1): 
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(This  is  a moment  per  unit  length  of  the  submarine  hull) 


The  next  two  contribution  integrals  (to  Eq.  9)  involve  the  fixed  vector  CQ 
which  was  defined  by  Eq  (5) : 
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The  results  of  Eq  (16)  and  (17)  represent  - as  was  the  case  of  Eq  (14)  - con- 
tributing moments  per  unit  length  of  the  submarine  hull. 
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Comparison  of  Complete  Analysis  to  an  Engineering  Approximation 

For  a submarine  of  the  "Trident"  class,  one  of  the  project  contractors 
intends  to  use  a heeling  moment  computation  that: 

(a)  Takes  into  account  possible  vertical  shear  of  the  current  im- 
pinging upon  the  hull  by  means  of  a quadratic  function  described 
in  our  Fq  HO). 

(b)  Computes  the  heeling  moment  of  the  submarine  by  a formula  that 
effectively  models  the  submarine  hull  as  a flat  plate  (perpen- 
dicular to  the  current),  of  area  equal  to  the  projected  area  of 
the  actual  submarine  hull. 

It  was  item  (b)  that  motivated  the  subject  research  note,  with  the  intuitive 
notion  that  the  above  approximation  (flat  plate)  is  overly  conservative  in 
the  sense  that  it  would  predict  too  large  a heeling  moment.  The  question 
remained:  by  how  much?  In  absence  of  current  shear,  this  question  can  be 
answered  by  a simple  parametric  analysis,  as  follows: 

Assumptions:  Both  analyses  use  the  same  C ^ (drag  coefficient)  for 

a plate  element  under  normal  incidence  of  flow. 

The  heeling  angle  is  (for  small  angles)  linearly  proportional  to  the 

heeling  moment  and  the  ratio  of  the  two  computed  moments  (per  unit 

length  of  submarine  hull)  will  be  the  ratio  of  the  resulting  angles. 

Under  this  assumption  both  analyses  are  performed  with  the  same  function 
G,  defined  by  Eq  (6).  Using  terminology  of  Fig.  1: 


"Flat  Plate"  Analysis 
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(z  + p ) Gdz 


(18) 


Eq  (18)  can  be  compared  to  the  slim  of  Eq  (14),  (16)  and  (17)  using  the 
following  geometrical  relations  (from  Fig.  1) : 
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Then,  the  "flat  plate"  analysis  yields: 
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When  the  same  relations  are  applied  to  Eq  (14)  thru  (17),  the  complete 
analysis  yields: 

"complete  analysis" 
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It  is  easy  to  see  that  for  small  a (refer  to  Fig.  1);  i.e. , for  a fully 
loaded  submarine,  the  ratio  of  Eq  (23)  to  (22)  is 


-1 


M 


co 


M. 


Co  F.P. 


1 + cos  a 


2 cos  3a 

3 + cos  a - -j-j 


2 

— 

J 

3 

y 

h n - 

0. 

(24) 


for  a -*■  0 
with  n - 0.02,  and  R/D  - 1 


Thus,  indeed,  the  approximate  engineering  analysis  is  conservative,  in  fact, 
possibly  too  conservative  because  it  seems  that  for  dockside  operations  on 
the  "Trident"  type  submarines  even  small  heeling  angles  would  be  on1y  re- 
luctantly tolerated. 


Illustrative  Quantification 

The  analysis  refered  to  on  p.  (the  "flat  plate"  approximation)  ends 
with  an  estimate  of  the  heeling  angle: 


®F.P.  = t(m 
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b(2240)  GM  J 


(25) 


where  Res  is  the  force  resultant  on  the  submarine  hull,  and  (p^  + pc)  its 
moment  arm  relative  to  C0.  A is  the  submarine's  displacement  in  "lohgtons", 
TjM  is  the  metacentric  height  in  (feet)  and  2240  converts  tons  into  pounds. 
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Actually,  in  our  terminology,  with  L the  length  of  the  submarine: 
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where  is  the  "specific  weight",  = p g,  of  the  water  and  the  simplifica- 
tion cancelling  L and  R is  possible  considering  the  submarine  hull  as  partial- 
ly submerged  cylinder.  Going  even  further,  using  Eq  (6)  and  assuming  a to 
be  small: 


F.P. 


tan 
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CD  U2  ^2_ 
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v g GM 


(27) 


With  CD  a 1.2,  U = 2 kt  = 3.37  ft/s,  Dq/R  b 1.1  and  GM  « 1 ft,  we  get 

®p.p . = 8,4°  while  a better  estimate,  using  results  of  our 


Eq  (24)  is  0 = 5.6° 


This  completes  the  analysis  by  J.  G.  Dworski.  Some  results,  computed 
from  actual  currents  in  Hood  Canal,  by  L.  H.  Larsen  follow  and  are  shown  in 
the  figures  hereafter. 
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FIGURES 


In  the  following  we  show  the  results  of  applying  EQ  (27)  to  the 
observations.  The  parameters  are  as  listed  below  EQ  (27).  Half-hour 
averaged  data  was  used.  Because  the  half-hourly  average  data  underestimates 
peak  currents  the  results  tend  to  underestimate  the  instantaneous  roll  as 
given  by  (27).  Also,  note  that  the  angle  is  almost  linear  with  respect  to 
the  argument  in  (27).  Thus,  for  example,  if  we  halve  GM  the  angles  are 
doubled . 

The  cumulative  plots  are  labeled  with  the  current  meter  number  and  the 
start  time  for  the  deployment  of  the  instrument.  In  all  cases  the  valid 
data  from  the  entire  deployment  is  used. 
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